Skip to content

Use bitset in VTKHDF read mesh, removing MPI_Allgatherv call. #3764

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Draft
wants to merge 6 commits into
base: main
Choose a base branch
from

Conversation

chrisrichardson
Copy link
Contributor

Rather than communicating all mesh cell types globally using MPI_Allgatherv use a bitset to do a set union on the cell type and degree. Using 8 cell types, and degrees 1-8, this fits into a 64-bit bitset.

@chrisrichardson
Copy link
Contributor Author

Addresses issue #3701

@garth-wells
Copy link
Member

I think a medium term fix would be that the mesh construction/factory function takes only the cell types on the calling process. Then this issue (largely) disappears.

Once the mesh topology is built, information on the cell polynomial degree can be communicated (at the same time that the geometry data is communicated.

@chrisrichardson chrisrichardson added this to the 0.11.0 milestone Jul 11, 2025
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants